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Abstract 

We consider the quantum-mechanical decay of a Schwarzschild-like black hole into 
almost-flat space and weak radiation at a very late time. That is, we are concerned 
with evaluating quantum amplitudes (not just probabilities) for transitions from initial 
to final states. In this quantum description, no information is lost because of the black 
hole. The Lagrangian is taken, in the first instance, to consist of the simplest locally 
supersymmetric generalization of Einstein gravity and a massless scalar field. The quantum 
amplitude to go from given initial to final bosonic data in a slightly complexified time- 
interval T = r exp(—i9) at infinity may be approximated by the form const. exp(— I), where 
/ is the (complex) Euclidean action of the classical solution filling in between the boundary 
data. Additionally, in a pure supergravity theory, the amplitude const. exp(-J) is exact. 
Suppose that Dirichlet boundary data for gravity and the scalar field are posed on an 
initial spacelike hypersurface extending to spatial infinity, just prior to collapse, and on 
a corresponding final spacelike surface, sufficiently far to the future of the initial surface 
to catch all the Hawking radiation. Only in an averaged sense will this radiation have an 
approximately spherically-symmetric distribution. If the time-interval T had been taken 
to be exactly real, then the resulting 'hyperbolic Dirichlet boundary-value problem' would, 
as is well known, not be well posed. Provided instead ('Euclidean strategy') that one takes 
T complex, as above (0 < 9 < n/2), one expects that the field equations become strongly 
elliptic, and that there exists a unique solution to the classical boundary-value problem. 
Within this context, by expanding the bosonic part of the action to quadratic order in 
perturbations about the classical solution, one obtains the quantum amplitude for weak- 
field final configurations, up to normalization. Such amplitudes are here calculated for 
weak final scalar fields. 

1. Introduction 

Since 1976, the most generally accepted option for the end-point of gravitational 
collapse to a black hole was that quantum coherence or information would be lost [1]. 
The gratifying change in opinion this past week [2] now makes it possible at last to begin 
publishing our late-1990's work on another option, namely, that information is not in fact 
lost, and that the end state is a combination of outgoing radiation states (see option 3 
below). The present work beg doctoral dissertation at Cambridge of one of the 

present authors [3] in 1997-2001 (the thesis was approved in summer 2002). This brief 
letter is an introductory sketch for numerous further papers contained in the thesis, which 
will provide detailed derivations of quantum amplitudes in black-hole evaporation. 
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For simplicity, we fix attention in this paper on non-rotating black holes. In the 
semi-classical theory, a black hole of initial mass Mq slowly radiates away its mass into 
predominantly massless particles, over a timescale ~ (Mq) 3 in Planck units [4,5]. Remark- 
ably, at late retarded times, the emission depends only on the late-time quasi-stationary 
state of the 'black hole' and not on the precise details of its formation. The frequency 
spectrum has precisely the thermal character expected, with the Hawking temperature at 
infinity (units restored) given by 



8nGM k B ' 

where ks denotes Boltzmann's constant. 

If one allows for the back-reaction of the emitted particles on the gravitational field, 
then, to a good approximation, one can model the black-hole evaporation with a sequence 
of quasi-stationary Schwarzschild solutions until the mass has nearly reached the Planck 
mass [4,5]. The Schwarzschild mass-loss rate balances the rate at which energy is radiated 
to infinity. 

For the (ill-understood) very late stages of this evaporation, there are a number of 
conceivable options [1,6-10]. Some of the main options are: 

1. The gravitational collapse could leave behind a naked singularity with negative 
mass (as there is no event horizon) which persists, thereby violating cosmic censorship; 

2. A stable Planck- mass object might remain, containing the negative-energy quanta 
and stellar matter, so that the joint quantum state of the remnant and the emitted quanta 
would be pure; 

3. The object might disappear completely, but all information about black-hole 
configurations and any locally conserved quantities might still escape to infinity. That is, 
particle correlations would be restored, so that the final state of the quantum field would 
be pure; 

4. The object might disappear completely, taking with it the information about the 
stellar matter, negative-energy quanta and any quantities not coupled to long-range fields, 
so that the final state would be mixed. The nature of the final radiation in this instance, 
and in option 3 above, is still to be uncovered. 

This paper focusses on options 3 and 4 and probabilities or quantum amplitudes asso- 
ciated with such outcomes. Below, we will describe how the total amplitude of (probable) 
weak-field final configurations can be computed from the imaginary part of the Lorentzian 
action, which is quadratic in the field perturbations in the linearised theory. The computa- 
tion of such amplitudes for late-time scalar perturbations is the main result to be reported 
in this paper. 

For simplicity, we take the bosonic matter to be described by a massless scalar field (p. 
The full Lagrangian might, in principle, be more complicated, containing also appropriate 
fermionic fields and (say) Maxwell or Yang-Mills fields. Indeed, although we shall not study 
fermionic fields in this paper, it will be simplest (see Sec. 2) to assume that the Lagrangian 
is that of [ ], consisting of supermatter - a complex massless scalar field and spin - 1/2 
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partner - coupled to N = 1 supergravity, consisting of Einstein gravity with gravitino 
partner. In a linearised approximation, we expand the massless scalar field and the metric 
around a background classical field configuration, in which <fi is taken to be real, describing 
a (bosonic) spherically symmetric 'matter' collapse, in which only the scalar field acts as 
a matter source for gravity. For both scalar and gravitational fields, 'Dirichlet boundary 
data' (see below) are specified on a (edgeless) spacelike hypersurface Ej at an early time 
- roughly the moment of collapse - and on a final late-time spacelike hypersurface E / , 
long after the black hole has evaporated. The location of E/ will be such as to register 
all the evaporated radiation. Both surfaces extend to spatial infinity. 

To fix one's physical intuition, it may be imagined that the energy density of the 
initial scalar configuration is extremely diffuse, such that almost all of the mass is dis- 
tributed over radii vastly greater than the 'Schwarzschild radius' 2M . Further, the initial 
time- variation of the scalar and the gravitational fields can be taken to be extremely slow. 

Away from the black hole, typical wavelengths of perturbation modes are small 
in comparison with the local radius of curvature of the space-time - the characteristic 
length scale over which the background space-time changes significantly. Prior to the 
black hole's final disappearance, however, typical wavelengths of perturbation modes which 
will contribute significantly to radiation in the late-time quantum treatment are compara- 
ble in magnitude to the black-hole radius. Correspondingly, one finds that the radiating 
Schwarzschild-like Vaidya metric [11-13] models, to a good approximation, the background 
space-time of an evaporating black hole at late times, outside the central strong-field re- 
gion. Such a semi-classical description will need to be joined onto a comparable treatment 
of the more non-linear central regions. 

Nevertheless, there is an intrinsically 'non-classical' character to this formulation, 
since wave-like boundary- value problems are well known to be badly posed [14]. Anal- 
ogously, the Uncertainty Principle forbids the simultaneous specification of a dynamical 
variable and its conjugate momentum - precisely the data specified in a Cauchy problem. 
The price that we pay for our boundary-value formulation may be a lack of existence 
and/or uniqueness in the solution of the boundary- value problem, in the case of a real 
time-interval T at infinity; this is intimately connected with the occurrence of singulari- 
ties. But experience of other simpler such 'hyperbolic boundary- value problems', such as 
the flat-space wave equation, leads one to expect that there will be a (complex) classical 
solution, typically unique, provided that we rotate the proper-time separation at infinity, 
T, infinitesimally into the complex, a la Feynman, and indeed that a real Riemannian 
solution will be reached when T is rotated by a phase of tt/2 to a Euclidean proper dis- 
tance r = %T. There is a substantial mathematical theory of such strongly elliptic Dirichlet 
boundary-value problems in the complex, with existence, uniqueness and regularity prop- 
erties corresponding to the more familiar real elliptic case [15]. 

1.1 Existence for the full Einstein— scalar system 

A variety of spherically symmetric boundary-value examples in the complexified ver- 
sion of general relativity, coupled to a massless scalar field, are being worked out, and will 
then be published. One would like to show that (three-dimensional) asymptotically flat 
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initial and final boundary data (Dirichlet data) for the spherically-symmetric gravitational 
and scalar fields can be specified, together with (in the first instance) a purely Euclidean 
'time-interval' |T|, such that there is a smooth classical Riemannian infilling solution of the 
coupled Einstein/massless scalar field equations, agreeing with the given Dirichlet bound- 
ary data. This requirement is perhaps not unreasonable, since the coupled Riemannian 
field equations are 'elliptic modulo gauge', and elliptic equations tend to smooth out bound- 
ary data. Further, that, as the time-interval-at-infinity T is rotated back in the complex, 
towards the Lorentzian region, by a phase less than tt/2, while holding fixed the Dirich- 
let data on the initial and final surfaces, then the full complex classical solution varies 
smoothly with the phase, as does the classical action. 

In particular, if the scalar-field part of the Dirichlet data is chosen to be sufficiently 
strong, then the corresponding Lorentzian-signature solution, reached by a phase rotation 
by 7r/2, will describe collapse to a black hole. Only as one approaches a phase of n/2 
would the classical solution become singular. Of course, this would be in accord with the 
singularity theorems of general Lorentzian-signature relativity [16]. This view of strongly 
elliptic versus singular Lorentzian-signature solutions also shows further the insight of 
Feynman's +ie prescription. 

Note that the mass M of the classical 'space-time' is a functional of the boundary 
data, including the possibly complex time-at-infinity T [17,18]. It has been found in the 
above examples that M is negative when T = —it, for r real, corresponding to a real 
Riemannian solution. The geometry is regular at the origin r = of spherical symmetry 
(and elsewhere) and the Riemannian solution lives on a portion of Euclidean space R 4 . One 
does not encounter the different topology of the Riemannian positive-mass Schwarzschild 
solution [19,20]. 

As described in [17,18,21], each of the intrinsic spatial metrics hijj and hijp on £/ 
and Hp will have an intrinsic Arnowitt-Deser-Misner mass Madm [22,23], measured by 
the usual 1/r fall-off rate of the spatial metric near spatial infinity. But (see Sec. 4. 4 of 
[18]) in general, for the classical infilling metric g^, the obvious surfaces of constant time 
will be emitted in such a way that the mass M of the 4-dimensional geometry differs from 
the '3-dimensional mass' Madm- R is this which allows for the possibility that the mass 
M of the classical geometry may depend on the time-separation T at infinity, as well as 
on the bounding 3-metrics hijj, hijF, as in the previous paragraph. 

Note also: it is essential that the 3-metrics prescribed on £/ and Hp should each have 
the same intrinsic ADM mass Madm- Otherwise [17,18], any classical infilling 'space-time' 
will have these 3-surfaces badly embedded near spatial infinity. 

In the more general non-spherical black-hole problem to be studied in the rest of this 
paper, the slight complexification of T induces an imaginary part in the total (Lorentzian) 
action, which leads (below) to the quantum amplitudes and probabilities needed to examine 
options 3 and 4 above. 

2. The quantum amplitude 

Consider first quantum field theory for gravity coupled to bosonic and fermionic 
fields. In particular, take the case in which non-trivial bosonic data are specified on the 
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boundary surfaces, while for simplicity the fermionic boundary data are taken to be zero. 
The 'Euclidean' quantum amplitude to go between the prescribed initial and final data is 
given formally by a Feynman path integral. Naively, disregarding the problem of infinities, 
one would expect the amplitude to have the asymptotic form 

Amp ~ (A + hA 1 + h 2 A 2 + ...) exp(-I B /h) . (2) 

Here Ib is the real classical 'Euclidean' action of a Riemannian solution of the coupled 
Einstein/bosonic-matter classical field equations. [For simplicity, we assume that there is a 
unique classical solution.] Further, Ib and the loop terms Ao, Ai, A2, ... depend in principle 
on the boundary data. For our case of matter coupled to Einstein gravity, Ib, Aq, A\, ... 
also obey differential constraints connected with the local coordinate invariance of the 
theory, and with any other local invariances such as gauge invariance [24]. 

In particular, when the theory is also invariant under local supersymmetry, and when 
the fermionic boundary data are now allowed to be non-zero, the above semi-classical 
expansion (2) becomes extremely simple. For example, for N = 1 supergravity, one has 
[18,25]: 

Amp ~ A exp(-I/h) . (3) 

This is not just a formal expression, as in Eq.(2), but an exact statement, freed of infinities. 
The 'one-loop factor' Aq is in fact a constant, while I denotes the action of the classical 
solution, which includes both bosonic and fermionic parts. Thus, in this case, the classi- 
cal action is all that is needed for the quantum computation. A corresponding situation 
arises with ultra-high-energy collisions, whether between black holes [26], in particle scat- 
tering [27], or in string theory [28]. For more complicated field theories possessing local 
supersymmetry, such as supergravity coupled to supermatter [29,30], quantum amplitudes 
may well still be meaningful since free of divergences. Their semi-classical expansion is, 
however, expected to have the more general form (2). 

The difficulty of assigning meaning to quantum amplitudes in field theories which 
contain Einstein gravity and other lower-spin fields, but which are not invariant under 
local supersymmetry, might also encourage one to turn to string theory. But it is again 
true in string theory that only for the superstring models, which are invariant under local 
supersymmetry transformations, does one have meaningful finite amplitudes [31,32]. 

In the asymptotically-flat, spatially- 1R 3 context here, the purely Riemannian case 
above corresponds to a time-separation at spatial infinity of the usual rotated form T = 
—it, where r is a positive imaginary-time separation. Consider further the locally super- 
symmetric case, but take only the bosonic data to be non-zero (real) on the initial and final 
surfaces. Suppose that there is a unique infilling Riemannian classical solution. Following 
the standard route, one can study the (now complex) bosonic amplitude (2) or (3) as T is 
rotated through angles 9 from 9 = n/2 to 9 = +e, with 

T = exp(-z#) r . (4) 

Provided that there continues to exist a (now complex) bosonic classical solution to the 
Dirichlet problem, as 9 > decreases from ix/2 towards zero, the expression (2) or (3) 
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should still continue to give the form of the (analytic) quantum amplitude. In particular, 
this would occur if strong ellipticity held for the coupled Einstein/bosonic-matter field 
equations. 

3. The classical action and amplitude for weak perturbations 

3.1 The approximate space-time metric 

The classical background bosonic fields, the metric g^ v and (real) scalar field <p, 
are each taken to have a 'large' time-dependent spherically symmetric part and a 'small' 
perturbative part, which can be expanded out in terms of sums over tensor, vector and 
scalar harmonics [33-39] . Each harmonic is weighted by a function of the time- and radial 
coordinates (£, r). The 'large' or 'background' Lorentzian space-time metric can be put in 
the form [40,41]: 

ds * = _ e b(t,r) dt 2 + e a(t,r) dj ,2 + r 2 ^2 + ^Q d ^ _ ^ 

The 'background' scalar field is denoted by <!>(£, r). The classical spherically-symmetric 
part of the Einstein and scalar field equations will be as given in [40] ('Lorentzian') or [41] 
(Riemannian) for an exactly spherically-symmetric set of fields, except for an additional 
effective energy- momentum contribution < T F J F >, resulting from local space-time aver- 
aging of the contribution to the Einstein equations of terms quadratic and higher in the 
combined perturbations [42,43]. It will be assumed that < T FFF > is defined in such a 
way as to be spherically symmetric. 

3.2 Scalar field — harmonic decomposition 

Now consider purely bosonic weak classical perturbations around this dynamic back- 
ground. For simplicity of exposition, consider a 'small' real scalar perturbation, denoted 
by <j)^\ where = $ + (p^\ Because of the approximate spherical symmetry, one may 
expand the perturbation (/>W in the form 

4>W(t,r,6,<p) = -J27- „E' _ Y £m (n) R im (t,r) . (6) 

Here, Yg m (O) denotes the (£, m) spherical harmonic of [44] . The boundary conditions on 
the radial functions Re m (t,r) as r^O follow from the regularity at r = 0. Recall that, 
in the strongly elliptic case, all fields are analytic in the interior of the large cylindrical 
boundary formed by the initial and final surfaces together with a surface at large r. 

The perturbed scalar wave equation at late times, 

VV^W =0, (7) 
with respect to the background geometry, leads to the (£, m) mode equation 

[ e ( b -^/ 2 d r } 2 R im - (d t ) 2 R irn - (^\[d t (a-b)} (d t R em )-V e (t,r)R em = 0, (8) 
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where 

b(t,r) 

V t (t,r) = G — 



(9) 



is real and positive in the 'Lorentzian' case. Here, m(t,r) is defined by exp[— a(t, r)] = 
1 — [2m(t, r)]/r . This potential Vg(t, r) generalises the well-known massless-scalar effective 
potential in the exact Schwarzschild space-time [23,39], which vanishes at the event horizon 
{r = 2M} and at infinity, and has a peak near {r = 3M}. 

There is, of course, an analogous harmonic decomposition of the weak gravitational- 
wave perturbations about the nearly-spherically-symmetric background [33-39]. 

3.3 The classical action 

After detailed calculation [3,45], one finds that the classical Lorentzian action for a 
classical bosonic solution to our model, including the quadratic-order contribution of weak 
perturbations, can be written as 



(10) 



Here, for simplicity, it is assumed that the only non-zero (prescribed) boundary perturba- 



tions are in the intrinsic spatial metric and in the real part of the scalar field on 
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the final surface Ef . The corresponding perturbations in the initial boundary data hf£ 
and 4>^ on Ej have been taken to be zero. Of course, one could easily include them also. 
Here the linearized perturbation of the momentum canonically conjugate to the 

intrinsic spatial gij, given in [23]. Also, 

tt^i, = n<°>" V^ 1 ' (11) 

is the momentum canonically conjugate to (p^K Here, the future-directed unit 

time-like normal vector to the three-surface (here E/), in the 'Lorentzian' case. Further, 
in Eq.(10), M is the '4-dimensional mass' of the 'space-time', as in Sec. 1.1. Typically, M 
will differ from the ADM mass Madm found from the intrinsic 3-metrics hiji, h^p. 

3.4 The quantum amplitude for scalar perturbations on the final surface 

On the initial surface Ej , we take the initial data such that there is a quasi-stationary 
Schwarzschild-like background. On or near the final surface E / , long after evaporation, 
one again expects a nearly- Schwarzschild background; the characteristic time-scale over 
which the background geometry varies is very much greater than the period of a typical 
wave. In such cases, the radial wave functions R£ m (t, r) can be decomposed harmonically 
with respect to t, and can be conveniently normalised (below). Hence, near E/ say, one 
can write 



r L — 'Cm 



/oo 
dk a U m Rk£m(t, r) Y tm (Q) , (12) 
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where the {cikim} are real quantities. Spherical symmetry implies that, on Ey , 

Rke m (t,r) = R k e(r) (13) 

is independent of m, where Ru{r) is a real function. Note [see below, Eq.(20)] that the 
allowed values of k are discrete, labelled by an integer n = 0,1,2,..., and are given by 
k = k n = n 7r/ 1 T | . 

We describe briefly the normalisation of the Rkt(r). With the notation f = kr, the 
boundary conditions of regularity as r^O imply 

Rkt(r) oc rjg(r) ~ const. r i+1 , (14) 

where jg(r) are (real) spherical Bessel functions [46]. For a background geometry which is 
nearly Schwarzschild at large radius, one also has 

Rki(r) ~ z u e lkr s + z u e~ lkr s (15) 

as r — > oo , where, for large r , 



r* s ~ r + 2Mlog^— J-lJ (16) 

is the Regge- Wheeler 'tortoise' coordinate [23,33]. Here the {zu} are dimensionless com- 
plex coefficients which are determined by the regularity above at r = 0, subject to suitable 
normalisation [3]. The Zkt are related to the Bogoliubov coefficients [47,48], thus making 
contact with the original formulation of black-hole evaporation [4,5]. The detailed relations 
are described in the fuller works [3,45]. 

Given the above information in a neighbourhood of the final surface E / , one can com- 
pute the classical scalar contribution to the Lorentzian action S c i ass [h^ , <f>^] of Eq.(10), in 
the case that the time-interval T = re at infinity is rotated slightly into the complex. 
This classical scalar contribution will in general be complex. The Euclidean action is then 
defined via 

-^class i S c \ ass . (17) 

In the limit that 6>^0 (while keeping 6 strictly positive), one finds that the (bosonic) 
quantum amplitude has the form (say, in the case of a locally supersymmetric Lagrangian): 

Amp ~ const. exp(- J c i ass ) (18) 

apart from possible loop corrections [Eq.(2)] which will only come into play at or above 
the Planck energy. This amplitude has a (computable) oscillating part in each mode, 
multiplied by a product of Gaussians 



in the coordinates {aum} of Eq.(12) for the perturbed final scalar data. Here, because of 
the large-but-finite time interval |T|, the eigen-frequencies are of the form 

k = k n = . (20) 

4. Comments and further work 

By using a field-based amplitude in quantum gravity and by rotating the time in- 
terval T at infinity slightly into the complex: T — > r e , one expects to have a strongly 
elliptic classical boundary-value problem. For weak wave fields of different spin at the final 
late-time boundary E / , after the black-hole evaporation, one can compute the quantum 
amplitudes for different final configurations as exp(— J c i ass ), up to an overall factor. An 
explicit example of this for the scalar-field perturbations of our model was given in Sec. 3. 
[Full details will be given in [45] and further papers.] Analogous 'Gaussian' results for the 
spin-1 Maxwell field and for spin-2 gravitational- wave perturbations are contained in [3] 
and will also be submitted; similarly for the fermionic case. 

The late-time radiation can be viewed either in a field or in a particle representation, 
and has the usual interpretation in terms of thermal Hawking radiation. 

Work is also under way on numerical solution of the complexified boundary-value 
problem for the coupled Einstein-scalar field equations, with real data posed on the initial 
and final three-surfaces, and a complex time-interval T specified at spatial infinity, as 
considered above. This requires generalisation of the original Riemannian Einstein/scalar 
numerical boundary- value work of [41]. This should help one to understand, for example, 
the nature of any burst of radiation which might originate in the late stages of complexified 
'collapse'. 

It would also be very interesting if similar methods could be applied in the case that 
the spatial bounding surfaces are not just two copies of M 3 . 
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